Strong correlations enhanced by charge-ordering in highly doped cobaltates 
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We present an explanation for the puzzling spectral and transport properties of layered cobaltates 
close to the band-insulator limit, which relies on the key effect of charge ordering. Blocking a 
significant fraction of the lattice sites deeply modifies the electronic structure in a way that is shown 
to be quantitatively consistent with photoemission experiments. It also makes the system highly 
sensitive to interactions (especially to intersite ones), hence accounting for the strong correlations 
effects observed in this regime, such as the high effective mass and quasiparticle scattering rate. 
These conclusions are supported by a theoretical study of an extended Hubbard model with a 
realistic band structure on an effective kagome lattice. 



The layered cobaltate metals are famous for their re- 
markable electronic properties [T], ranging from a large 
thermoelectric response [2] charge-ordering [3l E] and 
puzzling magnetic behavior, to superconductivity when 
intercalated with water [S]. Because of the universal- 
ity of the physical properties throughout the cobaltate 
family, it is believed that most of these properties have 
their roots in the electronic structure of the C0O2 layers. 
Those are formed of edge-sharing octahedra, with the Co 
ions forming a triangular lattice. Either alkali (-earth) 
metals (Na, Li, Ca etc., as in Na2;Co02) or more com- 
plex building blocks (e.g. rocksalt BiO or SrO planes in 
misfit cobaltates [6] ) separate the layers and serve as elec- 
tron donors. Varying the composition of the intercalated 
compounds controls the doping x € [0, 1], resulting in a 
nominal Co^''~^-'+ (Sd^"*"^) valence in the low-spin config- 
uration within the t2g manifold. The latter is split by the 
local trigonal symmetry into one aig and two e' states. 
Strong Coulomb interactions have been documented ff] 
to occur within this orbital subspace. 

Remarkably, most of the signatures of strong corre- 
lations in cobaltates, such as specific-heat enhancement 
[5] , band narrowing [71 IMTT] or strong electron-electron 
scattering [H] are observed in the high-doping region 
X € [1/2,1] close to the band-insulating limit {x = 1). 
Besides intriguing magnetic behavior, charge dispropor- 
tionation is a prominent feature in this part of the phase 
diagram. In Na3/4Co02 clear evidences of separation of 
Co-ions into nonmagnetic Co'^'^-ions and magnetic ions 
with a mixed valence are present [131 114) . Importantly, 
this charge segregation is not restricted to selected dop- 
ings, but the concentration of Co'^"'"-ions is experimen- 
tally verified to steadily increase with doping starting 
from X = 1/2 [T5]. Recently, the nuclear-magnetic- 
rcsonance (NMR) experiments of Alloul et al. |4j revealed 
for X — 2/3 a charge-ordering pattern of cobalt ions with 
valency close to Co^"*", with the remaining Co'^'^'''"'" ions 
forming an effective kagome lattice (EKL). The rather 
complex three-dimensional unit cell has common fea- 
tures revealed within a previous theoretical electronic- 



structure study [16j . In general, multiple sodium-ordered 
phases are found in highly-doped sodium cobaltates [TBl - 

In this Letter, we demonstrate that many puzzling 
physical properties of cobaltates at high doping may 
be understood within a picture of correlated itiner- 
ant electrons moving in a charge-ordered background. 
Charge-disproportionation patterns of Co^^-ions with 
well-localized electrons are formed, restricting the re- 
maining carrier dynamics to an effective lattice with an 
open structure. Such a system is shown to be highly sen- 
sitive to correlations, both because of a prominent van- 
Hove singularity close to x — 0.7 (owing to higher-order 
hopping terms) and due to the strong effect of intersite 
Coulomb interactions in line with the high susceptibility 
to charge-density wave (CDW) formation near x = 3/4. 

To substantiate this view we consider a particular re- 
alization of charge ordering, namely the EKL detected 
at a; = 2/3 for Naa;Co02 up to room temperature [3]. 
Whereas the charge-ordering at x = 1/2 seems unique 
to that stoichiometry, x = 2/3 also marks the onset of 
Curie-Weiss behaviour and enhanced spin fluctuations, 
leading eventually to in-plane ferromagnetic order for 
3/4 < X < 0.9. Furthermore, it is a commensurate 
doping for CDW instabilities on the triangular lattice 
(as studied e.g. in Ref. [T^ within the Gutzwiller ap- 
proach) . Importantly, a recent calculation of the charge 
susceptibility in sodium cobaltate for the Hubbard model 
shows tendencies to the kagome instability at large x even 
when disregarding the specific sodium distribution [2U] . 
Note however, that disordered donor configurations can 
by themselves enhance the correlation effects in the C0O2 
layers [21], also magneto-polaronic excitations were sug- 
gested to play a role in this context [22] ■ In what fol- 
lows, we address the spectral and transport properties of 
cobaltates stemming from a charge-order-induced EKL 
(see Fig. [T]), where the sodium and oxygen degrees of 
freedom are integrated out. The EKL doping y is re- 
lated to the physical doping x on the triangular lattice 
by y = (4a; — l)/3. The full t2g manifold is included in 




FIG. 1: (color online) Top left panel: Sketch of the EKL, with blocked sites depicted by small spheres and the shaded region 
representing the unit cell. Panels (a-e): Unfolded spectral functions y4(k, oj) (according to Eq. ([3|) projected onto aig orbitals 
(a-d) and onto both aig and the two e'g orbitals (e) for Na2/3Co02 {T — 115K). Top panels (a-b): results for {U — V = 0) on 
the triangular lattice (a), and on the kagome lattice (b), with the side-panels displaying the corresponding fc-integrated spectral 
function. Bottom panels (c-e): results obtained within CDMFT for (7 = 5 eV, V = (c) and for f/ = 5 eV, V = 0.3 eV (d), 
(e). The white dotted lines in (a-d) represent the QP peak and HE hump extracted from the ARPES data in Refs. [?1 [TUl In 
(e) the black dotted lines correspond to the locations of the energy-distribution-curve (EDC) maxima of the spectral intensity 
seen at high energies in the bare ARPES data of Refs|9l [10] (with black bars showing the approximate width of the feature) . 
The quasiparticle and HE hump dotted lines in (e) are the same as in the other panels. 



the modelling, however explicit two-particle interactions 
are only treated for aig, since it is known that the e' 
bands lie below the Fermi level Ef for the high-doping 
regime. Our realistic model on the EKL is defined by an 
extended Hubbard Hamiltonian H = Hq + Hi^t with 
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where the hopping matrix t^ is parametrized up to the 
third nearest-neighbor (NN) from the band structure of 
Naa;Co02 within the local-density approximation. Or- 
bital indices are a = 1 for aig and a = 2, 3 for e'gi, e^g- 
The parameter A fixes the aig-e' renormalized crystal- 
field splitting. We choose t/ = 5 eV in line with exper- 
imental and theoretical determinations [3 US] and var- 
ied the NN intersite Coulomb interaction V up to 0.3 
eV. Decreasing U down to 3.5 eV provided no significant 
variation of the results. 

The band structure of the noninteracting realistic 
model consists of nine bands inside the Brillouin zone 
(BZ) of the EKL. To connect to the results of angle- 
resolved photoemission (ARPES) experiments we repre- 
sent the computed fc-resolved spectral function in the BZ 
of the underlying triangular lattice. This can be achieved, 
following the method of Ref . |21l by " unfolding" the nine 



bands in the supercell kagome-lattice BZ (SBZ) to the 
three bands in the normal BZ (NBZ) of the triangular 
lattice. The states |Ka) defined within the SBZ are 
projected onto the states |ka) defined in the NBZ, so 
that the spectral function projected on orbital a reads: 

A^ (k, u) = -- ^(ka|K/) Im Gw (K, to) (KZ' |ka) , (3) 
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where K € SBZ and k G NBZ. 
lattice Green's function, Gn'(K,aj) = 

[{uj + ^)5w - Hq^iv (K) - S/;/ (K, cj)] " is written in 
the basis |/) = |«, a) , where i runs over the three atoms 
of the kagome-lattice unit cell. 

To examine the spectral properties, we project A(K, u) 
onto the aig and e' orbitals and unfold the EKL bands 
into the NBZ. The obtained spectral functions are com- 
pared to the ARPES data of Nicolaou et al. for x = 0.7 
O [TU]. Although those data were obtained for mis- 
fit cobaltates (BiO and BaO planes), at present stage 
the spectra obtained for sodium cobaltates are qualita- 
tively very similar (see e.g. [3S], [H]). The analysis of 
the polarization-dependence of the ARPES spectra per- 
formed by the authors of Ref. [HI [TU] led them to identify 
two oig-like features: a low-energy quasiparticle (QP) 
peak close to e-p, and a high-energy (HE) hump around 
-0.4 eV, both shown as dotted lines in Fig. [T] 

Already for U = V = Q the theoretical spectral func- 
tion of the EKL is found to be in much better agree- 



ment with experiment than that of the original triangular 
structure without charge order, as shown in Fig. [l^, b. 
The HE hump identified experimentally, and completely 
absent in the triangular-lattice model, matches quite well 
with the HE feature around -0.4 eV in the EKL model. 
This HE feature originates here from aig-e' hybridiza- 
tions (cf. Fig. flk where the full i2g projection is shown) 
induced by the loss of the original full triangular symme- 
try. It is thus essential to keep the e'g states in the model 
to reproduce the spectral structure in this energy range 
around the F-point. In contrast, in Ref. [^[TCTIthe HE fea- 
ture was interpreted as an incoherent excitation with aig 
character of unknown type. However, it should be noted 
that the full local triangular symmetry with equivalent 
Co sites together with a compensation between the spec- 
tra of the two e' orbitals were assumed in the analysis 
of Ref. [HI Uniin order to extract an aig signal. Since tri- 
angular symmetry is broken in the charge-ordered phase, 
refinements in the polarization-dependent analysis of the 
ARPES data may be required. 

The full interacting model is solved by dynamical 
mean-field theory (DMFT) and its cellular-cluster ex- 
tension (CDMFT) [211 117] ■ A continuous-time quantum 
Monte Carlo solver within the hybridization-expansion 
approach |28j is used, in the implementation by Par- 
collet and Ferrero [IS]- The cluster impurity model 
is defined in terms of a natural basis of molecular or- 
bitals via |di) = (-2|1) + |2) + |3))/y6, Ida) = 
(|2) - \3))/V2, \u) = (|1) -f |2) + |3))/y3, where \t) , 
i — 1,2,3 denote states corresponding to the basis atoms 
of the kagome-lattice unit cell. In this representation 
the impurity self-energy matrix is diagonal with elements 
(Sd, Srf, E„) = (Eo - El, Eo - El, Eo + 2Ei), where Eq, 
El are the on-site and intersite self-energies, respectively. 

In terms of spectral properties, the on-site [/-term 
alone (V — 0) already yields a significant renormaliza- 
tion of the QP band: its bandwidth is reduced by a 
factor Z ~ 0.4 (compare panels (b) and (c) in Figjl]). 
Spectral- weight transfer takes place to an incoherent up- 
per Hubbard band above ep (Fig. lit) . But the strongest 
effect of U actually manifests itself in a high value of the 
inverse QP lifetime F = — ZImE(iO"'"). Figure K^ displays 
the dimensionless ratio T/k^T vs. T evaluated in DMFT. 
Remarkably, this ratio is significantly larger than unity 
for almost all temperatures and doping studied. Thus, 
long-lived coherent QPs only exist at very low T, due to 
electron-electron scattering. Extrapolation of our numer- 
ical results suggest that the coherent regime {T/kBT < 1) 
would be reached only for T < 20K at a; = 2/3, and the 
true Fermi-liquid behaviour (T/kBT ex T) at an even 
lower temperature scale. This is qualitatively consistent 
with the large values of the resistivity reported in this 
doping regime (e.g. p ^ 300^r2cm for Nao.riCoOa at 
T — lOOK [IJ, as well as with the huge scattering rate 
at low T [T2| and the low value of the saturation temper- 
ature of the susceptibility [3D] . The value of F is largest 
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FIG. 2: (color online) Temperature-dependence of F/fcaT, 
with r the inverse QP lifetime for the kagome cobaltate model 
with on-site Coulomb interaction U {V = 0) for different x. 
Inset: F/fcaT for the triangular cobaltate model including dis- 
tant hoppings at X = 3/4 (green, circles) and for the kagome 
lattice with only NN hoppings at s = 2/3 (magenta, stars). 



for X ~ 0.69, which is approximately the doping at which 
the van-Hove singularity lies at the Fermi level (see side- 
panel of Fig. fTJD). Notably, only the EKL model with 
realistic hoppings shows this characteristics, while sole 
NN hoppings (having no flat dispersions near ep at large 
x for i < 0) lead to a low scattering rate (inset of Fig.[2|. 
The behavior of the system is strongly modified when 
the intersite interaction, y > 0, is taken into account. 
Figure [3] displays the QP weight Z„ from CDMFT for 
the 'symmetric' molecular orbital of the EKL as a func- 
tion of X for increasing V. For x = 2/3, V has only 
moderate effect, although it does reduce Z hy a factor 
of two, bringing it close to the experimentally observed 
mass enhancement (m*/mLDA — 4). Some improvement 
in the comparison to ARPES at this doping is also found 
(Fig. flH,e). In contrast, close to x = 3/4, the system 
is highly sensitive to V. Indeed, this doping level corre- 
sponds to the commensurate doping y = 2/3 on the EKL, 
at which a nonzero V can induce a CDW instability to- 
wards an insulating state. We find this to happen at a 
rather small value of 14 = 0.1 eV (inset of Fig. [3]). This 
renders it possible that a very strongly correlated metal- 
lic state (with a small QP weight renormalized by charge 
fluctuations) found close to 2; = 3/4, is best described as 
a weakly doped CDW insulator. It is also worth mention- 
ing that X = 3/4 is the doping for the onset of magnetic 
ordering, and that transport anomalies were observed at 
this doping |31) . 

Note that previous works already showed the impor- 
tance of nonlocal correlations [32] |33] . From a 'molecu- 
lar orbitals' viewpoint, the CDW transition is driven by 
orbital polarization: the intersite interaction drives the 
symmetric orbital \u) towards half-filling (hence suscep- 
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FIG. 3: (color online) Doping-dependence of Z for the sym- 
metric molecular orbital \u) of the kagome cobaltate model at 
T = 397^. Crosses: single-site calculations, thus Z^ = Zd = 
Z\ triangles, circles, boxes: cluster calculations. All lines are 
just guides to the eye. Inset: the local spectral function pro- 
jected onto the a^g orbital with V ?A, x = 3/4. 



tible to a MIT), while the two asymmetric orbitals |di,2) 
are completely filled (with in total 5 = 3(1 -fy) electrons 
in each EKL unit cell). 

In conclusion, we have shown that several of the in- 
triguing physical properties of layered cobaltates at large 
X can be explained within a picture of electrons moving in 
a charge-ordered background. In particular, good quanti- 
tative agreement with ARPES spectra has been obtained. 
Although we have considered a specific realization for 
the charge-ordered phase in the form of a kagome struc- 
ture, we expect similar effects for other ordering patterns 
in which a large fraction of Co ions localize charge and 
thus become blocked for electron hopping processes. This 
blocking makes the system highly sensitive to correlation 
effects, in particular to intersite Coulomb repulsion. The 
proximity of a van Hove singularity has also been shown 
to play a key role. Interesting issues to be addressed 
in future work are whether such charge blocking can ac- 
count for the large thermopower observed in this doping 
regime |34j and whether it can provide a mechanism for 
the apparent violation of the Luttinger sum rule [2] • 
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